Scissors mode of a rotating Bose-Einstein condensate 
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A scissors mode of a rotating Bose-Einstein condensate is investigated both theoretically and 
experimentally. The condensate is confined in an axi-symmetric harmonic trap, superimposed with 
a small rotating deformation. For angular velocities larger than u>±/^/2, where u>± is the radial 
trap frequency, the frequency of the scissors mode is predicted to vanish like the square root of the 
deformation, due to the tendency of the system to exhibit spontaneous rotational symmetry break- 
ing. Measurements of the frequency confirm the predictions of theory. Accompanying characteristic 
oscillations of the internal shape of the condensate are also calculated and observed experimentally. 
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Bose-Einstein condensates rotating at high angular ve- 
locity exhibit spontaneous breaking of rotational symme- 
try |IJ. This phenomenon, which is the consequence of 
two-body repulsive interactions, shows up in the occur- 
rence of considerable deformations of the trapped atomic 
cloud in the plane normal to the rotation axis. These 
configurations have been recently observed experimen- 
tally using a nearly axi-symmetric harmonic trap, with 
a small deformation of the trapping potential rotating 
around the z axis 

Under such conditions the collective modes of the sys- 
tem exhibit new interesting features. In particular one 
mode in the transverse plane (orthogonal to the rota- 
tion axis) corresponds to a shape-preserving oscillation 
of the atomic cloud with respect to the principal axes xy 
of the rotating trap. The frequency oj of this scissors- 
type motion is much smaller than the mean transverse 
trap frequency lo±_ . This surprisingly low value originates 
from the fact that the restoring force of the oscillation is 
proportional to the small trap deformation, while the mo- 
ment of inertia of the condensate remains finite due to its 
considerable deformation. This represents a major differ- 
ence with respect to the traditional scissors mode [[| ^] 
of a non-rotating condensate where both the restoring 
force and the moment of inertia vanish in the limit of an 
axi-symmetric trap. 

The purpose of this work is to provide both a theo- 
retic and experimental investigation of the problem. We 
first derive the relevant equations describing the scissors 
mode. For a rotation frequency fl larger than the critical 
value uj±/^/2 , we show that the scissors mode frequency 
uj vanishes like the square root of the trap deformation. 
We then report on the experimental observation of this 
scissors mode, and we present results which confirm the 
theoretic predictions with good accuracy. 

The atoms of mass m are confined in the transverse 
plane by the potential 

V ± (x,y)=m(L0 2 x x 2 +L0 2 y y 2 )/2, (1) 
where oj 2 y = uP[_(l ± e) . The {x, y) coordinates in the 



rotating frame are deduced from the coordinates in the 
laboratory frame by a rotation fit and the trap deforma- 
tion e is positive. A simple description of the collective 
oscillations of a rotating condensate is provided by the 
hydrodynamic equations of superfluids evaluated in the 
rotating frame: 
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V ■ [p(v - « A r)] = , (2) 



V + gp — mv ■ (ft At) = , (3) 



where p(r, t) is the spatial density, v(r, t) is the veloc- 
ity field in the laboratory frame, and fl = fl u z (u z unit 
vector along the z axis). The harmonic confining poten- 
tial V(r) — V±(x,y) + muj 2 z 2 /2 is time independent in 
the rotating frame. The parameter g characterizes the 
strength of the interatomic interactions and is related to 
the s-wave scattering length a by g — Ai:h 2 a/m. These 
equations are valid in the Thomas- Fermi limit, where the 
so-called quantum pressure term can be neglected in (|^) 
Q. We shall be interested here in vortex-free solutions 
for which V A v = . 

Stationary solutions of these equations can be obtained 
in the form Vo(r) = a'V(xy) for the velocity field and 
pa{r) = [p — m (uj 2 x 2 + uj 2 y 2 + u 2 z z 2 ) /2] jg for the den- 
sity profile. The density profile has the form of an in- 
verted parabola whose parameters (fi , oj x , ui y ) are deter- 
mined in a self-consistent way as a function of fl @. The 
value of the parameter a is determined by the solution(s) 
of the cubic equation a 3 + a(w^ — 2 fl 2 ) + euP^fl = . It 
is related to the deformation S = (y 2 — x 2 ) / (x 2 + y 2 ) of 
the atomic cloud by the simple expression a — ~fl5 . In 
the present work we will be interested in the so-called 
normal branch |Q, corresponding to the stationary solu- 
tions which can be obtained by an adiabatic increase of 
the angular velocity of the trap, starting from Q = (see 
Fig. Q). For fl > luj_/\/2 these solutions exhibit large 
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FIG. 1: Deformation S of the cloud for the steady-state solu- 
tions of Eqs. (|^) and (^) as a function of the angular velocity 
of the trap for a fixed value of the trap deformation e . The 
dotted line represents the solution for e = . 



values of the cloud deformation (5 ~ 1) even if the trap 
deformation e is much smaller than 1. 

An important class of collective oscillations can be 
derived on top of such stationary solutions by look- 
ing for time dependent solutions of the form 5p(r) = 
ciq + a x x 2 + a y y 2 + a z z 2 + a xy xy and cTv(r) = V(a 2; x 2 + 
a y y 2 + a z z 2 + a xy xy) where a, and on are time depen- 
dent parameters to be determined by solving Eqs. (|^) 
and (^). In the linear limit, one can look for solutions 
varying in time like e -1 "'. The collective frequencies ui 
are then found to obey the equation 
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with the coefficients Ci given in the appendix. In [jD it 
has been shown that these collective oscillations are dy- 
namically stable when evaluated on the normal branch. 
The stability conditions for higher multipole oscillations 
have been studied in 0, where it has been shown that 
the normal branch becomes dynamically unstable against 
the production of such excitations for f2 > 0.8 lu± when 
e -► 0. 

For the two external branches (see Fig. Q), one finds 
that the coefficient cq in (Q) vanishes linearly with e as 
e — > . In this limit, one of the eigenfrequencies of ([|) 
tends to 0, while the others remain finite. The frequency 
oj of the corresponding "soft" mode is: 



c _ ± \Qoj 2 SI(uj 2 _ 
ci 



3t4 +iuo 2 n 2 + 2f2 4 



(•5) 



The positive value for to 2 corresponds to the normal 
branch while the negative value, which is a signature 
of dynamical instability, corresponds to the opposite 
branch. For the normal branch, the frequency u> of the 
low energy mode vanishes like the square root of the trap 
deformation e , as announced in the introduction of this 



paper. It is worth noticing that this frequency does not 
depend on the value of u) z . By looking at the correspond- 
ing form of the solution for Sp and 5v , one can show that 
it corresponds to a scissors mode, i.e. to the oscillation 
of the angle 9 between the axes of the condensate and 
the principal axes x, y of the trap. The oscillation of 9 is 
accompanied by an oscillation of the condensate deforma- 
tion 5 around the equilibrium value So- This oscillation 
of 5 is dephased by tt/2 with respect to the one of the 
angle 9 . For small trap deformations we find 



o 2 y/2 n 2 - , 



(6) 



The fact that the frequency of the scissors mode van- 
ishes when e = (see Eq. (|j)) reflects the rotational 
invariance of the Hamiltonian. This behavior deeply dif- 
fers from the scissors mode in a non rotating condensate, 
whose frequency is given by u)\_\[2 and does not vanish 
in the limit of a symmetric trap || . The result (|5|) holds 
for a small amplitude A9 of the oscillatory motion of 9 . 
To study motions with a larger amplitude, one has to 
solve numerically the hydrodynamic equations (||) and 
(^) . We find that the frequency of the motion decreases 
as A9 increases, as for a simple gravitational pendulum. 

We now turn to the experimental observation of this 
scissors mode. We use a 87 Rb gas in a Ioffe-Pritchard 
magnetic trap, with frequencies oj x /2tt = 182 Hz, uj y ~ 
uj x , and oj z /2tt = 11.7 Hz. The cloud is pre-cooled using 
optical molasses to a temperature ~ 100 /iK. The gas is 
further cooled by radio-frequency evaporation to a tem- 
perature around 50 nK, corresponding to a quasi-pure 
condensate with 10 5 atoms. We denote by to the time at 
which the evaporation phase ends. For time t > to, the 
atomic cloud is stirred by a focused laser beam of wave- 
length 852 nm and waist wq — 20/im, whose position is 
controlled using acousto-optic modulators || . This laser 
beam creates a rotating optical-dipole potential which is 
harmonic over the extension of the cloud. The superposi- 
tion of this dipole potential and the magnetic one creates 
a transverse trapping potential identical to (|l|) . The trap 
deformation e is proportional to the laser intensity II 
and can be adjusted between and 4 %. 

The displacement of the center of the trap due to grav- 
ity and slight asymmetries in the trapping geometry pro- 
duce an additional static deformation which has been es- 
timated to be ~ 1% by measuring the splitting between 
the center of mass frequencies along the x or y axes. This 
static anisotropy plays a minor role in the present study 
and has been neglected in the analysis above. 

The frequency Cl(t) of the stirrer is first varied linearly 
during the time interval (to,ti), starting from Q (to) = 
uptoO(ti) = 2ttx139Hz (so that fi(ii) - 0.76w x ). The 
stirring frequency then stays constant in the time interval 
(ti,t2)- At time ti we switch off the magnetic trap and 
the laser stirrer, allow for a 25 ms free-fall, and image the 
absorption of a resonant laser beam propagating along z. 
We measure in this way the transverse density profile of 
the atom cloud, which we fit assuming a parabolic shape. 
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FIG. 2: Variation of the angle 8 (in degrees) and of the con- 
densate deformation 5 with the stirring time r = t% — 1\. The 
oscillations of 9 and 8 are a signature of the scissors mode. 
The oscillation frequency lu/2tt ~ 11.6 Hz is much smaller 
than the transverse trap frequency io±/2n, in agreement with 
Eq. ^. The solid lines are a fit to the data with a damped 
sine function. The dashed line in the lower graph shows the 
stationary value 5 = 8q. 



We extract from the fit the long and short diameters in 
the transverse plane (hence the deformation S(t)), and 
the orientation 9(t) of the condensate axes with respect 
to the rotating frame of the laser stirrer. 

The excitation of the scissors mode arises directly from 
the small non-adiabatic character of the condensate evo- 
lution as Cl(t) increases, during the time interval (io,ti)- 
At time t±, the state of the condensate slightly differs 
from the steady-state corresponding to Q(ti), Conse- 
quently the state of the condensate still evolves in the 
rotating frame in the time interval (£1,^2), even though 
the characteristics of the stirrer do not change anymore. 

We have plotted in Fig. || the angle 9 and the con- 
densate deformation 5 as a function of the time r = 
ti — t±. These data have been obtained for a ramp- 
ing time ti — to — 360 ms and a spoon anisotropy 
e = 0.017(6). The angle oscillates with a frequency 
oj/2tt = 11.6 (±0.1) Hz, in good agreement with the value 
11.4 (±2.1) Hz expected from (||). The initial amplitude 
A9 is 43 degrees, and the oscillation is damped with a 
time constant ~ 110 ms. The deformation 6 of the con- 
densate also exhibits a small oscillatory motion with an 
amplitude AS = 0.16 around the mean value Sq = 0.51. 
This motion has the same frequency as that of 9, is phase 
shifted by ~ 7r/2, as expected from (||) for the scissors 
mode, and it is damped with a similar time constant. 

To confirm that the oscillatory motion shown in Fig. ^ 
indeed corresponds to the scissors mode described above, 
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FIG. 3: Variation of the square of the scissors mode frequency 
ui/2n with the laser intensity II (arbitrary units). The error 
bars indicate the statistical spread of the results obtained for 
various runs of the experiment. The dotted line is a linear fit 
according to the prediction of Eq. 0. 



we have measured the variation of the frequency u) as a 
function of the laser intensity, which is itself proportional 
to the trap anisotropy e . The data are plotted in Fig. 0. 
They clearly show the expected dependance ui 2 oc e . Ac- 
tually, this scissors mode constitutes a very precise way 
to measure e , once the frequencies u>± and f2 are known 
with sufficient precision. 

We have also compared the measured amplitudes 
At9(t 2 ) and ASfe) with the results from a numerical in- 
tegration of the equations of motion for the coefficients 
di and cxi when f2 is ramped from up to its final 
value fi(ti). The calculated results reproduce the behav- 
ior found experimentally. The amplitude of the scissors 
mode decreases when the ramping time t± — to or the trap 
anisotropy e are increased. The physical reason for this 
variation is clear: the evolution of the condensate during 
the time interval (to,ti) is closer to adiabatic following, 
and the condensate is left at time t\ in a state closer to 
the stationary state expected for = fi(ii). However we 
could not reach a strict quantitative agreement between 
the calculated amplitudes and the measured ones (typi- 
cal deviation of 50%). We think that this is due to the 
damping of the scissors mode, present in the experiment 
(see Fig. ||) and neglected in our simple theoretical model. 
A proper description of this damping could be obtained 
using a formalism similar to 0, where the damping of 
the scissors mode in a static trap is investigated. 

Our final study concerns the behavior of the ratio 
AS/ A9 at t2, where A<5 and A9 are the oscillation ampli- 
tudes for the deformation and for the angle respectively. 
The results, plotted in Fig. [|, show that AS/ A9 varies lin- 
early with u}/2tt with a slope of 3.2 x 10 [degree. Hz] -1 
(i.e. A5/(ujA9) ~ 2.9 x 10~ 3 s/radian 2 )k This is in good 
agreement with the prediction of Eq. (0) which gives an 
expected slope of 4. Ox 10~ 4 [degree.Hz] -1 . The small de- 
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FIG. 4: Variation of the ratio AS/A8 as a function of u)/2n . 
The point with an horizontal error bar corresponds to the 
average of three results with very similar scissors frequencies 
ui. For small uo , corresponding to e — > 0, the oscillation of the 
cloud deformation is negligible and the condensate motion is 
similar to that of a rigid body. The dotted line is a linear fit 
to the data. 

viation between predicted and experimental slopes may 
be due to nonlinearities originating from the relatively 
large amplitude of the excitation. This linear dependance 
of AS/A8 on the scissors frequency confirms the idea that 
for almost symmetric traps, where e — > and thus lj — > , 
the motion of the condensate is shape-preserving, similar 
to that of a rigid body. 

To summarize, in this paper we have presented the the- 
oretic analysis and the experimental observation of a new 
mode of a trapped rotating condensate. The existence 
of this mode relies upon the breaking of the rotational 
symmetry of the condensate caused by atom interactions. 
The rotating condensate can have a sizeable deformation 
even in the limit of an arbitrarily small stirring poten- 
tial, and this gives unique properties to the mode of in- 
terest, such as a frequency much smaller than the trap 
frequency. Here we have been mostly interested in the 



small oscillations of the condensate around its rotating 
steady-state. A natural extension of this work consists 
in studying the non-linear regime of the system, where a 
chaotic dynamics can emerge. 
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Appendix 

The values of the coefficients c, entering into (||) are: 

c 3 = -c4{8(l + SI 2 ) + 3 A 2 } , 

c 2 = a4{4(5-2e 2 ) + 16^ 2 + 16f2 4 

-4[<S 2 (2S1 2 - 1) - 3dfk] 

+2A 2 (ll-a 2 + 12 2 )} , 
ci = -^i{16[(2fi 2 -l)(2^ 2 -l + e 2 ) 

-5 2 (2 fl 2 - 1 + 3 e 2 ) - &ne(5 - 4 - 3 2 e 2 ] 

+4 A 2 [(13 - 5 e 2 ) + 8 f) 2 + 12 fi 4 

~& 2 (2Cl 2 - 1) + 3afie]} , 
c Q = wl{40A 2 [(2fi 2 -l)(2O 2 -l + e 2 ) 

-a 2 (2 n 2 - 1 + 3 e 2 ) - 5fte(5 - 4 &) - 3 & 2 e 2 }} , 

where we have introduced the reduced quantities fl = 
Q/u>± , A = lu z /luj_ and a — a/uj± . 
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